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PREFACE 


Molecular interactions are the key for understanding the structure 
and properties of liquids. The physico-chemical behaviour and nature of 
intermolecular interactions in the liquid mixtures have been a subject of 
active interest, and fascinating progress in the study of intermolecular 
forces has been made in the last few years. Despite the large number of 
measurements on binary liquid mixtures, very limited studies have 
included pressure as an experimental variable, and in only very few cases 
these studies have been done at pressures other than atmospheric. 

Recently, study of the high pressure physics of liquids has become 
one of the most important fields in physical chemistry and chemical 
physics. But, due to scarcity of experimental data at elevated pressures, 
such studies could not be carried out for binary liquid mixtures. For the first 
time Takagi et al [J Chem Therm, 12, 277 (1980); 14 , 1167 (1982); 16 , 
1031 (1984); 17 , 1057 (1985); 20 , 809 (1988)] carried out accurate 
experimental measurements of ultrasonic velocity and related 
thermodynamic properties of some binary liquid mixtures at elevated 
pressures. In the present work, well-established Flory’s statistical theory [ 
Flory et al, J Am Chem Soc, 86, 3507, 3515, 3563 (1964); 87, 1833, 1838 
(1965) & J D Pandey et al, J Chem Soc Faraday Trans I, 84 , 1853 (1988); 
J Am Chem Soc, 104 , 3299 (1982); J Chem Research (s), 344 (2001); J 
Mol Liq, 100 / 2 , 153 (2002)] has been modified and developed for high 
pressure. This pressure-dependent Flory theory has been used for the 
computation of various thermodynamic properties of pure liquids and 
binary liquid mixtures at elevated pressures. 



Chapters 

Chapter-1 

Chapter-2 

Chapter-3 

Chapter-4 

Chapter-5 

Chapter-6 

Chapter-7 

Chapter-8 


CONTENTS 1 


Pages 

Introduction and literature survey. 1-14 

Theoretical formalism of thermodynamic 15-24 
properties of liquid mixtures on the basis 
of Flory's statistical theory - pressure 
dependent studies. 

Computation of excess volume, density 25 - 39 
and thermal expansivity at higher 
pressures using Flory’s statistical theory. 

Pressure dependent studies of internal 40 - 47 
pressure and pseudo-Gruneisen 
parameter for binary liquid mixtures. 

Estimation of pressure dependent thermo- 48 - 59 
dynamic properties (C p , C P E , C v , y) of 
binary liquid mixtures. 

Theoretical evaluation of ultrasonic veloc- 60 - 72 
ity, isothermal compressibility and 
isentropic compressibility of liquid mixtures 
at elevated pressures. 

Non-linearity parameter of pure liquids- 73 - 87 
pressure variation studies. 

Study of non-linearity parameter of binary 88 - 94 
liquid mixtures at elevated pressures. 





In the liquid state, thermodynamics and statistical mechanics can be 
used as a tool to further our knowledge about the interactions in liquid 
mixtures. Thermodynamics constitutes a very subtle study in physics, 
chemistry and engineering, and has innumerable applications. Statistical 
mechanics combines the laws of statistics with the laws of mechanics 
applicable to atomic or molecular particles. Thermodynamics presents 
broad relationship among the macroscopic properties but it is not 
concerned with the quantitative prediction of these properties. On the 
other hand, statistical mechanics seeks to establish relationship between 
macroscopic properties and intermolecular forces through partition 
functions and is very much concerned with the quantitative prediction of 
bulk properties. The statistical methods do not give the information of 
actual behaviour of the system but that of the maximum probability. The 
connection between statistical mechanics and thermodynamics is provided 
by entropy. 

There are two main reasons for the large amount of experimental 
and theoretical work on the properties of liquid mixtures' 1 ' 6,43 " 475 . The first is 
that they provide one way of studying the physical forces acting between 
two molecules of different species. The second reason is the appearance 
of new phenomenon which are not present in pure substances. 

Molecular interactions are the key for understanding the structure 
and properties of liquids. Nature of intermolecular interactions and 
thermodynamic behaviour of liquid mixtures have considerable 
significance in thermodynamics. Intermolecular forces on the whole are a 
very wide subject in the study of liquid mixtures. Various detailed theories 
and discussions in this regard have been given by Hirschfelder (1) , 
Margenau (2) , and Kihara (?) Rowlinson (4> , Hildebrand and co-workers (5) and 
Moelwyn-Hughes (6) . 

The subject of intermolecular forces is important from the view-point 
of its applications in many branches of science. Many attempts have been 
made to understand intermolecular attractions since the designation of 
“regular-solution theory” by Hildebrand' 75 . Prigogine and co-workers ' 8 ' 115 


developed the cell model, and also gave a theory based on the explicit 
combination of the cell method and the “corresponding state theorem”. 
Another statistical mechanical theory which has gained prominency is the 
“significant structure theory" developed by Eyring and his group (12 ' 16) . This 
is a purely thermodynamic theory, and yielded excellent results for most of 
the properties in the case of pure liquids and liquid mixtures. But this 
theory is very cumbersome and required a number of molecular 
parameters for its application. The hard sphere models proposed by 
various workers (17 ' 23) have also been used, in recent past, to estimate 
various thermodynamic properties of pure liquids and few binary liquid 
mixtures. Flory et al (24 ’ 27) used a simple partition function proposed by 
Eyring and Hirschfelder to get reduced equation of state which is based on 
the van der Waals potential energy model. 

Out of various statistical mechanical theories used in past, Flory’s 
statistical theory (24 ' 27) has gained much popularity due to its simplicity and 
validity. It has been applied to pure liquids, binary and multicomponent 
liquid systems of molecular liquids, metallic liquids, alloys and ionic liquids. 
In most cases, this theory has been found to be superior. The simplicity 
and precision in prediction of various thermodynamic properties makes the 
Flory theory a subject of special interest. Coefficient of thermal expansion 
(a) and isothermal compressibility (p T ) are the only two parameters of pure 
components which are needed for the computation of many 
thermodynamic properties. Moreover, it has no adjustable parameter. 
During recent years (28 ‘ 41,102 'i 05) , this theory has been widely used for 
estimating a number of properties of a variety of liquid mixtures. 

Despite the large number of measurements on binary liquid 
mixtures, few studies have included pressure as an experimental variable, 
and in only very few cases these studies have been done at pressures 
other than the atmospheric. It is widely known that the pressure has a 
great influence on packing and other effects in the mixture, so such 
measurements are suitable for testing the various theories of the liquid 
state. In the past few years, the study of the high pressure physics of 
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liquids have become one of the most important fields in physical chemistry 
and chemical physics. 

Excess thermodynamic properties provide information about the 
molecular Interactions and macroscopic behaviour of liquid mixtures. 
Excess thermodynamic functions are sensitively dependent not only on 
the differences in molecular forces but also on the difference in the size 
and shape of molecules. Use of excess function results from the fact that 
statistical mechanical theories of liquid mixtures naturally lead to the 
prediction of theoretical values of excess functions. 

There has been considerable interest < 24 - 27 ' 29 ’ 33 ) jn the theoretical and 
experimental investigation of the excess thermodynamic properties of 
liquid mixtures during the last few years. From the theoretical point of 
view, excess properties at elevated pressures are of particular interest, as 
they allow to test the applicability of molecular theories to predict the 
thermodynamic behaviour of liquids and their mixtures at ambient 
conditions. 

The study of excess compressibilities and excess volume provides 
important information on the intermolecular forces determining the 
properties of liquid mixtures. Excess molar volume reflects not only 
interactions but also difference between free volumes of the components. 
It has been established that the sign and magnitude of V E give a good 
estimate of the strength of the unlike interactions in a binary mixture (42) . 

Ultrasonic velocity has been found to be the most powerful tool for 
the structural and interaction studies of liquids and liquid mixtures. 
Ultrasonic propagation parameters and thermodynamic parameters are 
very important in understanding the physico-chemical behaviour of the 
liquid state. Measurement of sound velocity provides a convenient method 
for determining those properties of liquids which are not easily accessible 
by direct experimentation. 

Measurements of ultrasonic velocity have been adequately 
employed in understanding the molecular interactions in pure liquids and 
liquid mixtures. The propagation of ultrasonic (sound) wave through a 


medium is essentially a thermodynamic property, and is related to the 
physico-chemical properties of the system. Recent developments have 
made it possible to use ultrasonic waves in medicine, engineering, 
agriculture and chemical industry. Ultrasonic studies find extensive 
applications in almost every field of science. 

Most of the physical, thermodynamic and acoustic properties of the 
liquid mixtures have been theoretically (47,48) derived from those of the 
components, by assuming that the dissimilar molecules are non- 
interacting. This case exists in mixtures where both the components are 
non-polar. In the binary mixture, where one component is polar and other 
is non-polar, interactions have been observed (49) . Even if the interaction 
between polar and non-polar molecules is weak, there will be considerable 
change in the molecular environment. In such a case, the physical, 
thermodynamic and acoustic properties are likely to be affected by 
intermolecular hydrogen bonding, dipole-dipole, dipole-induced dipole and 
charge transfer interaction between the unlike molecules due to the 
variation in the relative concentration of the components 149,601 . 

Successful attempts have been made in the recent past (51 ' 56) on 
theoretical evaluation of ultrasonic velocity and its correlation with other 
thermodynamic properties in binary liquid mixtures using statistical and 
semi-empirical theories. 

Sound velocity depends on the density of the medium. The density 
of liquid mixture and related volumetric properties are required for many 
applicative aspects, as well as for theoretical calculations. Particularly, 
their knowledge is indispensable when the classical and thermodynamic 
approach is to be used to relate and to compute the equilibrium properties 
of mixtures. From this point of view, they are of strategically importance for 
both the test of existing theories and the development of new models for 
mixture behaviour. 

Thermal expansion coefficient is a very important thermodynamic 
property and knowledge of this property helps in calculating various 
thermodynamic parameters. Recently (57) an empirical relation has been 


proposed for computation of thermal expansivity which uses only the 
values of ultrasonic velocity and density at given temperature. 

isentropic compressibility Is a highly useful property in 
understanding the nature of molecular interactions in liquid mixtures 
because it depends upon the volume, and volume is very sensitive to the 
intermolecular interactions in liquid mixtures. In recent years, there has 
been an increased interest in the determination of either isentropic 
compressibility or isothermal compressibility. The role that isentropic 
compressibility plays in the thermodynamic characterization of binary 
liquid mixtures is also a subject study in recent years (58) . The isentropic 
compressibility may be used to separate the contribution of pressure and 
temperature in the change of thermodynamic properties. 

The more fundamental thermodynamic quantity is the isothermal 
rather than the isentropic compressibility and is fairly sensitive to a change 
of pressure. It is the key parameter in molecular thermodynamics of fluid 
phase equilibria. 

Isothermal compressibility has been widely evaluated by many 
workers* 19 ' 21 * for pure liquids but very limited attempts* 1 7,59 * 61) have been 
made to evaluate isothermal compressibility for liquid mixtures. 

One of the most important properties of the liquid state is the 
internal pressure which is the resultant of the forces of attraction and 
repulsion between the molecules of a liquid. The internal pressure of 
liquids is highly useful in understanding molecular interactions, internal 
structure, clustering phenomenon and ordered structure. Attempts have 
been made by Barton (62) , Rosseinsky (63) , Berkowitz et al* 64) Hildebrand et 
al (65) and Suryanarayana (66) to demonstrate the significance of internal 
pressure as a fundamental property of the liquid state and its correlation 
with other thermodynamic properties. Theoretical deduction of internal 
pressure of binary liquid mixtures was carried out by Subrahmanyam et al 
* 67 ' 69) using Flory’s statistical theory. 


Another quantity of interest is the Griineisen parameter, which was 
originally defined for solids but has recently been shown to be useful in 
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studying liquids. The concept of Griineisen parameter has been extended 
to liquid state for investigating the internal structure, clustering 
phenomenon and the remaining quasi-crystalline lattice nature. The 
Griineisen parameter is a dimensionless measure of change in entropy 
with volume or thermal pressure. Several workers have evaluated internal 
pressure 167 ' 71 * and Griineisen parameter 172,73 * in binary liquid mixtures. But, 
as far as our knowledge is concerned, no work has, so far, been done on 
the pressure dependent properties of binary mixtures on the basis of Flory 
theory. 

Heat capacity measurements (74,75) have been found to be useful in 
the investigation of change of liquid structure taking place during mixing. 
Heat capacity is unique amongst thermodynamic quantities as an indicator 
to changes of structure (i.e. order) or non-randomness in pure liquids and 
mixtures during the mixing process. The order-destruction and order- 
creation processes in mixtures can also be studied through determination 
of the corresponding excess quantity. Heat capacity measurements were 
carried out by Patterson et al (76,77) and Ahluwalia et al (78,79) during recent 
years. Flory’s statistical theory has been utilized to estimate heat 
capacities in various binary mixtures but in only very few cases these 
estimation have been done at pressures other than atmospheric. 

It is well known that the finite-amplitude waves travelling through a 
fluid will give rise to a number of non-linear phenomenon such as wave 
distortion, harmonic generation and sound saturation. It has become of 
much interest to predict the extent of various kinds of non-linear effects, 
which are due to deviation from linearity in ordinary acoustics 180 ' 84 *. The 
acoustician who studies finite-amplitude sound propagation in liquids is 
naturally interested in the non-linearity parameter, B/A, of liquids. The non- 
linearity parameter is a basic parameter for determining the degree of 
waveform distortion. It plays a significant role in non-linear acoustics and 
its determination is of increasing interest in a number of areas ranging 
from underwater acoustics to medical science. From the knowledge of the 
parameter of non-linearity, one can gain information about some physical 



properties (85,86) of the liquids such as internal pressure, intermolecular 
spacing, acoustic scattering and structural behaviour etc. 

Study of non-linearity parameter of liquid mixture is important \ 
in view of the information it may yield on the interaction in liquid mixtures. 

B/A values of the liquids have been interpreted as the quantity 
representing the magnitude of the hardness of liquid (87) which may be 
considered to be true for the liquid mixtures as well. 

Beyer et al* 80 " 82 ) were the first who presented the theoretical 
formalism of non-linearity parameter and determined it for a number of 
organic liquids using sound velocity and other thermodynamic data. During 
recent years* 88 ' 92,96 ' 1 ° 0) , a number of experimental and theoretical studies 
have been performed on the non-linearity parameter of liquids and liquid 
mixtures. Also, measurement of this parameter for biological jj 

media* 93 " 95) including living tissues is one of the more recent interests as it | 

can provide important information for the ultrasonic application in 
biological research related to diagnosis and therapy. Although detailed I 

studies of the non-linearity parameter of liquids and liquid mixtures have | 

been reported, to the best of our knowledge, computation of B/A at 
elevated pressures has not been done so far. 1 

Thermophysical and thermoacoustical parameters have been | 

successfully used to study the intermolecular interactions in liquids. These 1 

properties have been expressed in terms of the volume expansion, which 
is found as the controlling factor for their dependence on temperature and 
volume for liquids. Beyer’s non-linearity parameter (B/A)* 81, 82) , Rao’s 
acoustical parameter (K)* 101) and the Grtineisen parameter (r)* 72,73) are all 
simply related to each other. All the thermoacoustical properties are 
obtained with the help of thermal expansion coefficient whose dependence 
on temperature and pressure may conveniently be utilized to obtain them. 

Thus, it is evident that the study of the properties of liquid' state 
consists of a wide range of problems and to discuss all of them will 
certainly not be exhaustive but practically beyond the scope of the present 
thesis. During the past few years, the study of the high pressure physics of 


liquids has become one of the most important fields in physical chemistry 
and chemical physics. In the present thesis, some important 
thermophysical and thermoacoustical parameters of pure liquids and liquid 
mixtures have been evaluated using Flory’s statistical theory at elevated 
pressures. As far as our knowledge is concerned, no work has, so far, 
been done on the pressure dependent properties of binary liquid mixtures 
on the basis of Flory theory. This is our first attempt to apply Flory theory 
for evaluating properties of binary liquid mixtures at elevated pressures. 

The present investigative work entitled “Statistical Mechanics and 
Thermodynamics of Liquids and Liquid Mixtures” is compiled in 
forthcoming eight chapters. The entire thesis work encompasses 
theoretical predictions and interpretations of thermodynamic and other 
properties of liquid mixtures. 

Chapter-1 of the present thesis is introductory. It is devoted for the 
objective, scope and outline of the present work. An upto-date literature 
survey alongwith the background of the work is presented in this chapter. 
At the end of this chapter brief outline of the work presented in the 
subsequent chapters are given. 

Recently, extensive work on Flory theory has been done in 
substantial amount at atmospheric pressure. In chapter-2 of the thesis, 
Flory theory has been developed for high pressure to evaluate various 
thermodynamic properties of liquids and liquid mixtures. 

In chpater-3 of this thesis, the values of excess volume, density and 
thermal expansion coefficient of two binary liquid mixtures e. g. toluene + 
o-xylene and toluene + aniline, have been computed at elevated pressures 
and at temperature of 303.15 K using Flory’s statistical theory by the 
method described in the previous chapter. 

The aim of chapter-4 is to evaluate internal pressure and pseudo- 
Grwneisen parameter of two binary liquid mixtures using pressure- 
dependent Flory’s statistical theory at temperature 303.15 K. The binary 
mixtures under consideration are toluene + o-xylene and toluene + aniline. 


Chapter-5 of the thesis deals with the computation of heat capacity 
at constant pressure, heat capacity at constant volume and heat 
capacities ratio of the above mentioned binary liquid mixtures using Flory 
theory at temperature 303.15 K and at pressures ranging from 0.1 MPa to 
160 MPa. Excess heat capacities of the liquid mixtures have also been 
reported. 

Chapter-6 comprises the theoretical predictions of ultrasonic 
velocity, isothermal and isentropic compresibilities of the two binary 
mixtures viz. toluene + o-xylene and toluene + aniline at elevated 
pressures and at temperature 303.15 K for four different mole fractions 
using Flory theory. 

Chapter-7 deals with the computation of acoustic non-linearity 
parameter, B/A, of pure liquids e.g. toluene, o-xylene and aniline at 
elevated pressures and at temperature of 303.15 K using different 
approaches, The different methods used are Ballou’s empirical relation, 
Hartmann’s relation, Tong-Dong method using Schaaffs equation for 
sound velocity and Sharma’s method. A comparative study of B/A values 
obtained from different methods has also been carried out. 

Chapter-8 consists of the theoretical evaluation of non-linearity 
parameter of binary liquid mixtures. The mixture under present 
investigation are toluene + o-xylene and toluene + aniline. The methods 
used and the conditions are the same as employed in the previous 
chapter. 
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Chapter-2 


Theoretical formalism of thermodynamic properties 
of liquid mixtures on the basis of Flory’s statistical 
theory - pressure dependent studies 



Introduction : 


There has been an increasing interest in the study of 
thermodynamic functions as they are useful in understanding the nature of 
intermolecular interactions between component molecules of a mixture. 
Among all the theories, known so far, the Flory’s statistical theory (1 ' 6) has 
been the most widely applied theory of liquid mixtures. It has been 
employed to describe the thermodynamic behavior of binary and 
multicomponent systems by various workers {7 " 18) . But, it appears from 
literature survey that most of the studies have been limited to zero 
pressure* 7 " 18 *, and no attempt has been made, so far, to evaluate the 
thermodynamic properties of pure and binary liquid mixtures at elevated 
pressures. 

For the first time Takagi et al (20 ' 24) carried out experimental 
measurements of ultrasonic velocity and related thermodynamic properties 
of some binary liquid mixtures at elevated pressures. 

The present work deals with the theoretical evaluation of 
various thermodynamic properties of pure liquids and binary liquid mixtures 
at elevated pressures. Here the reduced equation of state (4) due to Flory is 
solved to get expressions for the thermal expansion coefficient (a) and 
isothermal compressibility (p T ) under the restriction of pressure 
considerations . Further, pressure - dependent expressions are derived for 
the reduced volume (?) and reduced pressure (?) which are used to get 
the different thermodynamic properties of pure components and binary 
mixtures namely, toluene + o-xylene and toluene + aniline. 

Theory : 

According to Flory’s statistical theory (4) , based on van der 
Waals potential energy model, the reduced equation of state is given by 

PV F l/3 1 


T 


( 2 . 1 ) 
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where the reduced parameters v ,f and Pare given by 

V -V IV* 

T=T/T’ 

and, P = p/p > 

(2.2) 

Rearranging eq. (2.1) , we get, 

x_(V ,n -l )(PV 2 + l) 

y 4/3 

(2.3) 

or, 


T=PV+V~' i -PV 213 -V~ An 

(2.4) 


Differentiating eq.(2.4) with respect to f keeping P constant, we have 


In), 


4~ - 

b- T/ 

3 


+ -F - 7 
3 


• (2.5) 


Again, differentiating eq. (2.4) with respect to p keeping f constant and 
rearranging, we obtain, 


'€) fj 

dpU 


2 +| PV~ 


3-p1 = k2/3^1/3_ 1 ) ....(2.6) 


The coefficient of thermal expansion (a) and isothermal compressibility 
(j3 r ) may be expressed in terms of reduced variables as (1) , 


(dv s 

\ _ ? \ 


tar. 

y .. 
■;^ii 

^1 

m 


■ (2.7) 


and , 



•( 2 . 8 ) 


Now using eqs. (2.5) and (2.6) in eqs. (2.7) and (2.8) respectively, and 
solving to eliminate f in each case, we get 


(aT)~ ] 


2 PV 2 1 

(PV 2 +1) 3(F 1/3 -l) 


(2.9) 


and 




i 


- 1+2 


PV Z J [3(K -1 

Solving eqs. (2.9) and (2.10) together, and rearranging we obtain 
ccT 


-(2.10) 


y_\ aT T 
[ + 3(l - 2f} r P + aT )_ 


( 2 . 11 ) 


Using this value of v is eq. (2.9), we get 
1 


js_ PtP 
(ptr-fiS) 


h 3(l - 2 J3 r P + aT) 


•( 2 . 12 ) 


Thus, with the knowledge of thermal expansion coefficient and isothermal 
compressibility, we can calculate reduced volume (V) and reduced 

pressure (P) for pure components at elevated pressures with the help of 
eqs. (2. 11) and (2.12). Equations (2.11) and (2.12) are the pressure- 
dependent expressions for the reduced volume and reduced pressure. 
Reduced temperature (f) can be calculated from eq. (2.3) using the values 
of V and P . Characteristic parameters can be obtained with the help of 
eq. (2.2). 

Binary Liquid Mixtures : 

Segment fraction (y/)and site fraction (<9)of liquid mixtures 
have been computed by the method suggested earlier (8,9,11 ' 18) according to 
which 


where V i is the characteristic volume of the i th pure component and 


0i = - 


Wt + 2 y k 

k*i 



-( 2 . 14 ) 


thus for binary mixtures, the segment fractions and site fractions will be 


¥2 ~ 


X/2 

xiK+x 2 V 2 * 


and i// i = 1 — if / 2 


( 2 . 15 ) 


and, 


¥2 


V 2 

vi+.n * 


*V/3 


and 9 { =\-0 2 


The interaction parameter x u is defined as (8 



( 2 . 16 ) 


( 2 . 17 ) 


where p' and pj are characteristic pressures of pure 
components i and j. 

For binary liquid mixtures, there is only one Flory’s interaction parameter 
term X 12 which is expressed as, 



The characteristic pressure P* for the mixture is calculated by 
the following equations (8,9,11,18) 


p = 




£ Vi QjXjj 


-(2.19) 


thus, for binary mixtures 

P + ^2 )~Vl^42)] 


The characteristic temperature r* of the mixture is given by 18 ' 9 ' 11 ' 18 * 


-( 2 . 20 ) 


Vfit 


-( 2 . 21 ) 


I 

4 T, 


For binary mixtures 7** will be 
T* -- 


vA | ViP* 
71* T‘ 


and the reduced temperature f of the mixture is given by 

~ T 
T 

The ideal reduced volume Pis given by 
F°=Ix,F f 
For binary mixture 

V° — +x 2 V 2 

The ideal reduced temperature f°is defined with the help of eq. (2.3) 

--(2.26) 


-( 2 . 22 ) 

(2.23) 

(2.24) 

- (2.25) 




0^3 


The excess reduced volume, v * , is defined as (2,5) 


dV ■ (f-f°) 
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using eq. (2.5), with negligible approximation, we get 


F=|p_p- J _ |/T°"*+| f°' K ] \t-T°) --(2.27) 

The reduced volume, V , for the mixture is given by 

V = V°+V E (2.28) 


with this value of reduced volume V of the mixture, various parameters for 
the liquid mixture can be calculated. 

Excess Volume (V E ) : 

The excess volume of the mixture is calculated with the help 
of equation (i0) . 

— (2.29) 


Thermal Expansion Coefficient (a) : 

Thermal expansion coefficient of the mixture is calculated 
with the help of pressure - dependent expression given by 


1 2 ^ 2 l 


(2.30) 


Isothermal Compressibility (J3 r )\ 

The pressure-dependent expression for the isothermal 
compressibility is given as 



(2.31) 


Internal Pressure (/>) : 

Internal pressure is related to the thermal expansion 
coefficient and isothermal compressibility through the relation, 


P, = 



(2.32) 


Heat Capacity at Constant Pressure (C,,) : 

Heat capacity at constant pressure of the mixture is defined 


as, 

C P = C p (idl)+ Cp 


(2.33) 


where C p (idl.) is the ideal heat capacity of the mixture and is defined as, 


C P (idL)= Zx/C/3 (2.34) 

where C Pi is the heat capacity of i th pure component 

The excess heat capacity, C/; , is defined by Khanwalkar et al. (12) according 

to Flory theory as, 



For binary mixtures, eq.(2.35) reduces to 


P*v' 

1 

1 ' 

T* 




Adiabatic Compressibility (j3 s ) : 

From the well-known thermodynamic relation 


Pt-P* 


a TV 
C„ 


— (2.35) 


(2.36) 


(2.37) 


we get, 


Ps 



■(2.38) 


with the help of this equation, we can calculate the value of adiabatic 
compressibility for the liquid mixtures knowing a, p T and C P . 


Heat Capacities Ratio (y ) : 

The well-known thermodynamic relation for the heat capacities 
ratio is given as, 


r = 


Pt 

Ps 


(2.39) 


Heat Capacity at Constant Volume (C K ) : 

The heat capacity at constant volume of the mixture is 
calculated by the relation 


Cy 


£jl 

r 


(2.40) 


Pseudo - Grwneisen Parameter (r) : 


Pseudo-Grwneisen parameter of the mixture can be 
expressed as (19) 



(2.41) 


Density (p ) : 

The molar volume (V) of the mixture is related with the 
characteristic volume (F‘)and reduced volume (V) of the mixture as, 

V = V*V ..... (2.42) 

where 

V * = X x,V* 


(2.43) 
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and Fis obtained from eq. (2.28). 

We know that density is defined as 

P = ~f (2.44) 

where M m is the molecular weight of the mixture given by 

M m = I x t M i — (2.45) 

where Mj is the molecular weight of the i th pure component. 

Thus, with the help of eq. (2.44), we can calculate the theoretical values of 
density for the mixture. 


Ultrasonic Velocity(w) : 

Ultrasonic velocity of the mixture can be calculated with the 
help of equation , 


u = 



— (2.46) 
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Introduction: 

The deviations from ideal behaviour can best be expressed in 
terms of excess thermodynamic functions* 1 *, which give more quantitative 
idea about the nature of molecular interactions. The difference between 
the thermodynamic function of mixing for a real system and the 
corresponding value for an ideal system at the same temperature and 
pressure is called thermodynamic excess function. It is denoted by 
superscript E. This quantity represents the excess (positive or negative of 
a given thermodynamic property of the solution over that in the ideal 
solution, e.g. 

X 6 =X-X i( ii 

Excess thermodynamic properties 12 ’ provide information about 
the molecular interactions and macroscopic behaviour of liquid mixtures, 
and can be used to test, and to improve thermodynamic models for 
calculating and predicting fluid phase equilibria. Excess thermodynamic 
functions are sensitively dependent not only on the difference in 
intermolecular forces but also on the difference in size and shape of 
molecules. 

During the last few years (3) , the interest in experimental and 
theoretical work on thermodynamic properties of liquid mixtures at 
elevated pressures has greatly increased. From the theoretical point of 
view, excess properties at elevated pressures are of particular interest, as 
they allow to test the applicability of molecular theories to predict the 
thermodynamic behaviour of liquids and their mixtures at these conditions. 

Values of excess molar volume provide a very sensitive test 
for the molecular solution theory. It has been observed that the V E of 
mixing of two liquids can arise from any one or more of the following 
factors: 

(I) difference in size and shape of the component molecules, 

(II) structural changes, such as changes in the correlation of 
molecular orientation, 
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(III) differences in the intermolecular interaction energy between like 
and unlike molecules, and 

(iv) formation of new chemical species through hydrogen bonding 
or electron donar acceptor interactions. Normally, volume changes occur 
because of the combination of these volume factors. 

Excess molar volume, V E , reflects not only interactions but 
also difference between free volumes of the components. The excess 
volume can be divided into two terms (5) ; V s f ze , associated with the size 
difference of molecules and the other Vj„ t , associated with the interaction 
forces between molecules. Thus, for a liquid mixture, we have 

V E = v? + v.e 

v v size v int 

The term V s f ze is negative if the radii of both components are different, 
and is zero if the radii are equal. The values of are observed to be 
negative or positive. This term is regarded as the excess volume which 
occurs from an interaction between molecules. 

It has been established that the sign and magnitude of V E give 
a good estimate of the strength of the unlike interactions in a binary 
mixture (6) . Large positive excess volume are taken as indicative of weak 
intermolecular interactions where as large negative values of excess 
volume are usually found when the interactions are strong and 
intermolecular complexes are believed to form. Extensive work on excess 
volume of liquid mixtures has already been carried out and still in progress 
by various workers (7 ' 12) . 

A number of thermodynamic properties can be evaluated 
when the sound velocity data are combined with density, pressure and 
temperature. Ultrasonic velocity, in conjunction , with density 
measurements, permits the direct estimation of adiabatic compressibility 
(J3 S ) and other related parameters which are not easily accessible by any 
other means. 
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Thermal expansion coefficient (a) is a very important 
thermodynamic property and knowledge of this property helps in 
calculating various thermodynamic parameters. Recently <13) , an empirical 
relation has been proposed for computation of thermal expansivity which 
uses only the values of ultrasonic velocity and density at given 
temperature. 

In the present work, the values of excess volume, density and 
thermal expansion coefficient of binary liquid mixtures have been 
computed at elevated pressures using Flory’s statistical theory by the 
method described in chapter two. For the prediction of thermal expansion 
coefficient of liquid mixtures, the reduced and characteristic parameters 
are calculated with the help of thermal expansion coefficient and 
isothermal compressibility of pure components. Molecular weights of the 
pure components are taken from the literature <4) . Here, we have used 
recently proposed relations 03,15 * for computation of thermal expansion 
coefficient and isothermal compressibility of pure components, 


a - 


75.6 *10" 3 

T m u \n p m 


(3.1) 


Pt ~ 


17.1* IQ" 4 

r 4 / V 4/ V 


(3.2) 


The other parameters of pure components are calculated as follows: 

Isentropic compressibility (/? s ): Isentropic Compressibility have been 
calculated using the well-known thermodynamic relation, 

p s = 4 - (33) 

up 

where the experimental values of ultrasonic velocity and density of 
pure components are taken from the paper of Takagi et al (14) 
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Heat Capacity at Constant Pressure^) :lt is given as, 

a 2 TV 


C P =- 


Pt~Ps 


-(3-4) 


Heat Capacities Ratio (x): The well-known thermodynamic relation for 
the heat capacities ratio is given as, 


/ = @JL 
Ps 


-(3.5) 


Heat Capacity at Constant Volume ( C v ): The heat capacity at constant 
volume is calculated by the relation, 

Cy=^- (3.6) 

r 

Pseudo- Gr^neisen Parameter(r) : Pseudo- Grwneisen parameter can 
be expressed as, 

r =^ (3J) 

Internal Pressure (/)): Internal pressure is related to the thermal 
expansion coefficient and isothermal compressibility through the relation, 

Pi=~-P (3.8) 

Pt 

All these parameters of pure liquid components are reported in 
table 3.1. 

Theory: 

Equations (2.29), (2.30) and (2.44) have been used to obtain 
the values of excess volume, thermal expansion coefficient and density 
respectively of the binary liquid mixtures under consideration. 
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Results and Discussion: 

The values of excess volume, density and thermal expansivity 
of two binary mixtures have been evaluated with the help of Flory’s 
statistical theory at elevated pressures, and at constant temperature of 
303.15 K. The two binary liquid mixtures under the present consideration 
are: toluene + o-xylene and toluene + aniline. 

The density and ultrasonic velocity of pure liquid components 
and binary mixtures are taken from the literature 04 ’. Thermal expansivity 
and isothermal compressibility of pure liquids are evaluated from relations 
(3.1) and (3.2). The desired characteristic and reduced parameters have 
been calculated using the methods discussed in chapter two. Ail the 
parameters of pure liquid components are taken from table 3.1. Tables 3.2 
and 3.3 enlist the computed values of excess volumes of binary mixtures 
toluene + o-xylene and toluene + aniline respectively. Tables 3.4 and 3.5 
enlist the experimental and computed values of density alongwith the 
percentage deviations at different compositions of toluene + o-xylene and 
toluene + aniline respectively. Tables 3.6 and 3.7 contain the experimental 
and computed values of thermal expansion coefficient alongwith the 
percentage deviations at different concentrations for the respective 
mixtures. 

A close perusal of table 3.2 shows that the computed values 
of excess volume for toluene + o-xylene are, generally, positive in 
magnitude with very few exceptions at higher pressure side, indicating 
weak interactions between component liquids. But in the case of toluene + 
aniline (table 3.3) the values of excess volumes are found to be negative, 
which show strong interaction between component molecules of the 
mixture. 

From tables 3.4 and 3.5, it is clear that computed values of 
density show the same trend as experimental one for both liquid mixtures 
under the present investigation. For both the liquid mixtures, the tables 
show that as the mole fraction of toluene in the mixture decreases, the 
density of the mixtures increases for the entire range of pressure. It is also 
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clear that for each mole fraction, the value of density increases as the 
pressure increases. The percentage deviations are found to be less than 
one for toluene + o-xylene, and for toluene + aniline it is found to be 
somewhat greater but maximum value is below ten percent. Tables 3.6 
and 3.7 reveal that the computed values of thermal expansion coefficient 
for both the mixtures toluene + o-xylene and toluene + aniline show same 
trend as experimental one. For both the liquid mixtures, it is clear that as 
the mole fraction of toluene in the mixture decreases, thermal expansivity 
decreases over the entire range of pressure. Also, the values of thermal 
expansion coefficient decrease as the pressure increases at all mole 
fractions with a very few exceptions in the case of toluene + aniline. 
Percentage deviations are less than one for toluene + o-xylene and less 
than three percent with some exceptions for toluene + aniline system. 



Table 3.1 
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Calculated values of excess volume of Calculated values of excess volume of 
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Table 3.4 
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Introduction: 

Internal pressure has been found to be a very important 
parameter in the theory of liquid state. In liquid state, cohesive forces are 
of prime importance. Different liquid systems differ essentially in their 
relative degree of cohesion and this cohesion creates a pressure directed 
inward of magnitude between 10 3 to 10 4 atm. Thus, the internal pressure 
which is the resultant of forces of attraction and repulsion between the 
constituents in the liquid medium provides the estimate of pressure 
directed inward by cohesion. The internal pressure of liquids is highly 
useful in understanding molecular interactions, internal structure and the 
clustering phenomenon. 

Pioneer attempts by Hildebrand et al (1) , Barton (2) , 
Rosseinsky (3) , Berkowitz et al (4) and Suryanarayana (5) demonstrate the 
significance of internal pressure as a fundamental property of the liquid 
state and its correlation with other thermodynamic properties. 

Internal pressure has been a subject of active interest during 
recent past. Successful attempts have been made by several 
• investigators (6 ' 10) to determine the internal pressure of binary liquid 
mixtures. Theoretical deduction of internal pressure of binary liquid 
mixtures was carried out by Subrahmanyam et al (6 " 8) using Flory’s 
statistical theory. 

In the theory of lattice dynamics, the extent of anharmonicity of 
the lattice is characterized by a parameter called Grwneisen parameter 
denoted by T. This dimensionless parameter is basically a diagnostic 
parameter and has been found to be a very useful tool for investigating 
the thermodynamic and other anharmonic behaviour of the solid 
crystalline lattices (11 " 15) . 

For the first time Knopoff and Shapiro (16) extended the study 
of this parameter to liquids also. It has been found that Gri/neisen 
parameter is an equally suitable parameter to investigate the internal 
structure, clustering phenomenon and other quasi-crystalline properties of 
liquids. 
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For liquids, the pseudo-counterpart of the term r is defined as, 

r = au 2 /C P (4.1) 

Where u is the velocity of sound, a the coefficient of thermal 
expansion and C P is the heat capacity at constant pressure. 

Using the relation, u 2 = — — 

P S P 


, Cp f$T 
and, y = — L - = £ - L ~, we get 

c r Ps 
r aV _ r- 1 
Pt C V otT 

where V is the molar volume. 


.(4.2) 


Several workers (17 ‘ 18) have evaluated Grwneisen parameter in 
binary liquid mixtures and this parameter has been used to study the 
intermolecular interactions between the molecules. 

In the present work, we have computed internal pressure and 
pseudo- Grwneisen parameter of two binary mixtures namely toluene + o- 
xylene and toluene + aniline using pressure-dependent Flory’s statistical 
theory. As far as our knowledge is concerned, no work has, so far, been 
done on the pressure dependent properties of binary mixtures on the basis 
of Flory theory. This is our first attempt to apply Flory theory for evaluating 
internal pressure and pseudo-Griineisen parameter of binary liquid 
mixtures at elevated pressures. 

Theory: 

Equations (2.32) and (2.41) have been utilized to obtain the 
internal pressure and pseudo- Grwneisen parameter respectively of liquid 
mixtures on the basis of Flory’s statistical theory as discussed in chapter - 
2 of this thesis. 
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Results and Discussion: 

The two binary liquid mixtures under consideration for the 
present investigation are: toluene + o-xylene and toluene + aniline. 
Internal pressure and pseudo-Grwneisen parameter of the two above 
mentioned liquid mixtures have been evaluated at elevated pressures 
ranging from 0.1 MPa to 160 MPa and at temperature of 303.1 5 K for four 
different mole fractions. Flory’s statistical theory has been employed to 
compute the internal pressure and pseudo-Grwneisen parameter of the 
binary mixtures at the above mentioned conditions. Experimental data of 
pure components and binary mixtures have been taken from the 
literature (19) . Wherever experimental values are not available, these are 
obtained with the help of well-known empirical relations recently (20 ' 22) 
proposed using the experimental values of density and ultrasonic velocity 
available in the literature (19) . Theoretical results, thus obtained, have been 
compared with the experimental findings. The experimental and calculated 
values of internal pressure, along with percentage deviations at different 
mole fractions for toluene + o-xylene and toluene + aniline, are reported in 
tables 4.1 and 4.2 respectively. Similarly, the experimental and theoretical 
values of pseudo-Griineisen parameter, along with percentage deviations, 
for the aforementioned mixtures are reported in tables- 4.3 and 4.4 
respectively. 

A careful study of tables 4.1 and 4.2 show that the calculated 
values of internal pressure show the same trend as observed 
experimentally. In all the cases, the values of internal pressure increase 
with the increase in pressure. The values of internal pressure increase 
with decrease in mole fraction of toluene in both the liquid mixtures at all 
pressures. The percentage deviations between computed and observed 
values of internal pressure are below one percent for toluene + o-xylene 
and less than two percent for toluene + aniline. 

But, the calculated values of pseudo-Gri/neisen parameter 
show opposite trend than that observed experimentally for both the liquid 
mixtures. 


Table 4.1 
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Experimental and calculated values of internal pressure of x-toluene + (1-x) aniline at 303.15 K 
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Introduction: 

Heat capacity measurements have been found (1,2) to be useful 
in the investigation of change of liquid structure taking place during mixing. 
Heat capacity is a second - order thermodynamic quantity and it seems to 
be unique amongst thermodynamic quantities as an indicator to changes 
of structure (i.e. order) or non-randomness in pure liquids and mixtures* 3 * 
during the mixing process. Thus, order-destruction and order -creation 
processes in mixtures can also be studied through determination of the 
corresponding excess quantity. The negative excess heat capacity (Cp) 
corresponds to order-destruction process whereas for order-creation 
process, the excess heat capacity (Cp ) is positive. 

Heat capacity measurements were carried out by Patterson et 
al* 4,5) and Ahluwalia et al (6,7) during recent years. Flory’s statistical theory 
has been utilized to estimate C P and Cp in various binary liquids. Despite 
the large number of measurements on binary liquid mixtures, few studies 
have included pressure as an experimental variable, and in only very few 
cases these studies have been done at pressures other than, atmospheric. 
In the present work, we have applied Flory’s statistical theory at elevated 
pressures. 

Theory: 

In the present chapter C P , C P , C v and y of binary liquid 
mixtures have been calculated with the help of Flory’s statistical theory at 
elevated pressures through the eqs. (2.33), (2.36), (2.40) and (2.39) 
respectively. The reduced and characteristic parameters, which are used 
in computation, have been deduced by the procedure mentioned in 
chapter-2. Parameters for pure components are obtained from table 3.1. 
Thermal expansion coefficient and isothermal compressibility of pure 
components are obtained from empirical relations* 9-11 * recently proposed 
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and isentropic compressibility is obtained from the thermodynamic 
relation, 

P s - (u 2 P )~ ] . where the values of ultrasonic velocity and density are 

taken from the literature (8) . Heat capacities at constant pressure are 
obtained from the relation, 

c a2w 
Pt ~Ps 


Results and Discussion: 


The two binary liquid mixtures taken for the present 
investigation are : toluene + o-xylene and toluene + aniline. Heat capacity 
at constant pressure, heat capacity at constant volume and heat 
capacities ratio have been evaluated using Flory’s statistical theory at 
elevated pressures at 303.15 K. The calculated values for both the 
mixtures are compared with the experimental ones, and percentage 
deviations were calculated. Excess heat capacities of the liquid mixtures 
have also been reported. Since for the above mentioned liquid mixtures, 
the experimental values of C P , C v and y were not available, we have 


adopted an indirect method to get these values. The well-known 
thermodynamic relation is, 


Pt ~Ps 


a 2 TV 

C P 


Rearranging this equation, we get 

C P =J^ 

Pt ~ Ps 


so, we can write 


C P(exp.) ~ 


g 2 (exp.)TF (exp ) 

PTfexp,) ~ Ps(exp.) 


,(5.1) 


•(5.2) 

.(5.3) 


where, isentropic compressibility 
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^(exp.) 2 

U (exp.)/?( exp .) 


..(5.4) 


The experimental values of ultrasonic velocity and density are taken from 
the literature {8) . Thermal expansion coefficient and . isothermal 
compressibility are obtained from the empirical relations^' 1 1} using 
experimental values of ultrasonic velocity and density at temperature T as, 
75.6 *10" 3 


* (exp) “ 


....(5.5) 


. „ . 17.1 *1(T 4 

and, p (exp.) = -— 1 — 

r 4/ V(exp.)p 4/3 (exp.) 


..(5.6) 


Experimental heat capacities ratio is obtained from 

0T (exp .) 

^(exp .) = ~ 7T 

(exp .) 

and experimental heat capacity at constant volume as, 


C F(exp.) 


^(exp.) 

y (exp.) 


....(5.7) 


.(5.8) 


Tables 5.1 and 5.2 enlist the calculated and experimental 
values of heat capacity at constant pressure along with the percentage 
deviations for the systems toluene + o-xylene and toluene + aniline 
respectively. Tables 5.1(a) and 5.2(a) record the values of excess heat 
capacity for mixtures toluene + o-xylene and toluene + aniline respectively. 
The experimental and calculated values of heat capacity at constant 
volume along with percentage deviations are reported in tables 5.3 and 
5.4 respectively for toluene + o-xylene and toluene + aniline. In tables 5.5 
and 5.6, the experimental and calculated values of heat capacities ratio 
along with percentage deviations at four mole fractions are reported 
respectively for toluene + o-xylene and toluene + aniline. 



A close perusal of tables 5.1, 5.2, 5.3, 5.4, 5.5, and 5.6 
reveals that calculated values of C P , C y , and 7 for both the mixtures 

show the same trend as experimental one. The values of C P are, in 
general, negative for both liquid mixtures. A close observation of tables 5.1 
and 5.2 show that C P increases as pressure increases at every mole 
fraction and at a constant pressure, C P increases as the mole fraction of 
toluene in the mixtures decreases. For toluene + o-xylene, the maximum 
deviation is nearly five percent but in most cases it is less than one 
percent. But for toluene + aniline, the maximum percentage deviation is 
about ten. In this case percentage deviation increases as mole fraction of 
toluene in the mixtures decreases. 

From tables 5.3 and 5.4, it is evident that variation in C v with 


pressure and also with mole fraction show, in general, the same trend as 
C P i.e. Cy increases when pressure is increased or the composition of 
toluene in the mixtures is decreased. Percentage deviations for C v are 
greater than that for C P , Percentage deviation increases as pressure 
increases for toluene + o-xylene but for toluene + aniline, this variation is 
not so straight forward but some what complex. A careful study of 
tables 5.5 and 5.6 reveals that value of y decrease as pressure increases 
at every mole fraction for both the mixtures. But at a constant pressure 
when mole fraction is changed, an interesting result is seen for both the 
mixtures. The value of /decreases as content of toluene decreases in the 
mixtures upto a certain pressure but it increases as toluene decreases 
after this certain pressure. This pressure is 30 MPa for toluene + o-xylene 
and 70 MPa for toluene + aniline. For both the mixture^ percentage 
deviation increases with pressure and maxijj^^^SMagie elation is 
upto twenty for both the mixtures. 

2 . 0 . 





Table 5.1 
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Table 5.3 
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Table 5.4 
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Table 5.5 
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Table 5.6 


58 



59 


References: 

1. LA Ball, D Patterson, M Costas & M Caceres, J Chem Soc Faraday 
Trans I, 84, 3991 (1988). 

2. LA Ball, M Costas, P Paquet, D Patterson & M E St. Victor, Pure 
Appl Chem, 61, 1075 (1989). 

3. M Costas & D Patterson, Thermochim Acta, 120, 161 (1987). 

4. A Dominguez, G Tardajos, E Aicart, S Perez-Casas, L M Trejo, M 
Costas, D Patterson & H Tra, J Chem Soc Faraday Trans, 89, 89 
(1993). 

5. H Van & D Patterson & J Soln Chem, 1 1 , 793 (1 982). 

6. J C Ahluwalia, R Bhat & N Kishore, J Chem Soc Faraday Trans I, 
84, 2651 (1988). 

7. J C Ahluwalia & N Kishore, J Soln Chem, 19, 51 (1990). 

8. T Takagi & H Teranishi, J Chem Therm, 17, 1057 (1985). 

9. J D Pandey, G P Dubey, N Tripathi & A K Singh, J Int Acad Phy S.ci, 
1,117 (1997). 

10. J D Pandey & R Verma, Chem Phys, 270, 429 (2001 ). 

11. J D Pandey, J Chhabra, R Dey, V Sanguri & R Verma, Pramana J 
Phys, 55(3), 433 (2000). 



Theoretical evaluation of ultrasonic velocity, isothermal 
compressibility and isentropic compressibility of liquid 
mixtures at elevated pressures. 



Introduction: 

Ultrasonic velocity has been a subject of active interest 
during the recent past as it plays an important role in the investigation of 
structure and intermolecular interactions between components of liquid 
mixtures. Ultrasonic propagation parameters (ultrasonic velocity, 
absorption, adiabatic compressibility etc.) and thermodynamic parameters 
(isothermal compressibility, thermal expansion coefficient, heat capacities, 
heat capacities ratio etc.) are very important in understanding the physico- 
chemical behavior of the liquid state. In addition, ultrasonic velocity has 
been proved to be an important property to test the validity of liquid state 
models. 

There are certain thermodynamic properties of liquid which 
are not easily accessible by direct experimentation. Measurement of 
sound velocity provides a convenient method for determining these 
thermodynamic properties of liquids. Literature survey (1 " 8) shows that 
ultrasonic study of liquid mixtures is highly useful in understanding the 
nature of molecular interactions present in the mixtures. Various 
successful attempts have been made in the recent past (9 ' 14> on theoretical 
evaluation of ultrasonic velocity and its correlation with other 
thermodynamic properties in binary liquid mixtures using statistical and 
semi-empirical theories. But such studies have been limited to binary liquid 
mixtures at zero pressures. Due to the lack of experimental data on 
ultrasonic velocity at elevated pressures, such studies could not be 
extended at elevated pressures. For the first time, Takagi et al (15 " 18) have 
carried out accurate experimental measurements of ultrasonic velocity and 
related thermodynamic properties of few binary liquid mixtures at elevated 
pressures over a wide range of temperatures. Although of late (19 " 27) , some 
work on ultrasonic speed and thermodynamic parameters has been 
carried out, however, to our knowledge, comparatively limited work is 
available at elevated pressures. In the present work, an attempt has been 
made to evaluate theoretically ultrasonic velocity at elevated pressures 
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using Flory’s theory. The theoretical method employed has already been 
detailed out earlier. 

Isentropic compressibility (p s ) has been widely used to study 
the molecular interactions through its excess value (39 " 41> . On the other 
hand, it can also be used to deduce other useful thermodynamic 
properties e. g. isothermal compressibility (p T ), heat capacities ratio (y), 
internal pressure (Pj) etc. These properties are accessible from isentropic 
compressibility (p s ) provided the values of thermal expansion coefficient 
(a) and heat capacity at constant pressure (C P ) are known. Isentropic 
compressibility is a very complex quantity. The more fundamental 
thermodynamic quantity is the isothermal rather than the isentropic 


compressibility. The isothermal compressibility, ’ is a 

directly measurable thermodynamic quantity. It is the key parameter in 
molecular thermodynamics of fluid phase equilibria. The isothermal 
compressibility of a liquid or liquid mixture can be obtained directly from 
pressure-volume measurements, or indirectly from the isentropic 
compressibility (p 8 ) according to the relation 


where a is the thermal expansion coefficient, T the temperature, V 
the volume and C P is the heat capacity at constant pressure. 

Isothermal compressibility has been widely evaluated by 
many workers (28 ' 32) for pure liquids but very limited attempts (33 ' 36) have 
been made to evaluate isothermal compressibility for liquid mixtures. In 
the present work, Flory’s statistical theory has been employed for the 
theoretical evaluation of ultrasonic velocity, isothermal and isentropic 
compressibilities of binary mixtures at elevated pressures. 
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Theory: 

Equations (2.46), (2.31) and (2.38) have been used to obtain 
ultrasonic velocity, isothermal compressibility and isentropic 

compressibility respectively of the liquid mixtures under the present 
consideration with the help of Flory’s statistical theory at elevated 
pressures. The reduced and characteristic parameters, which are used in 
computation, have been deduced by the method discussed in chapter -2. 

Results and Discussion: 

The two binary liquid mixtures taken for the present 

investigation are: toluene + o-xylene and toluene + aniline. The values of 
ultrasonic velocity, isothermal compressibility and isentropic 

compressibility of the binary mixtures have been evaluated with the help of 
Flory’s statistical theory at elevated pressures ranging from 0.1 MPa to 
160 MPa and at constant temperature of 303.15 K for four different mole 
fractions. Parameters of the pure components are obtained from table-3.1. 
Thermal expansion coefficient and isothermal compressibility of pure 
components are obtained from the empirical relations* 37,365 where the 
values of ultrasonic velocity and density are taken from the literature 065 . 

The calculated values of the properties for both the mixtures 
are compared with the experimental ones and percentage deviations were 
calculated. Experimental values of ultrasonic velocity for both the mixtures 
are taken from the literature 065 . Since the experimental values of 
isothermal and isentropic compressibilities were not available, we have 
adopted an indirect method to get these values. Using well-known 
thermodynamic relation, 

A = 4- -(6.1) 

UP 

we can write Ar e x p.) =: ”1 

U (exp.)/? (exp.) 


--( 6 . 2 ) 
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The experimental values of isothermal compressibility are obtained 
using empirical relation (37,38) . 


Pt\ 


(exp.) ' 


17.1 * 1 0~ 4 

T 4/9 M 2 (exp.)>0 4/3 (exp.) 


-(6.3) 


The experimental values of ultrasonic velocity and density of liquid 
mixtures, employed in the above eqs. (6.2) and (6.3), are taken from the 
literature (16) . 

The experimental and calculated values of ultrasonic velocity 
along with the percentage deviations at different mole fractions for toluene 
+ o-xylene and toluene + aniline are reported in tables 6.1 and 6.2 
respectively. Tables 6.3 and 6.4 contain the experimental and calculated 
values of isothermal compressibility with percentage deviations at different 
mole fractions for toluene + o-xylene and toluene + aniline respectively. 
Similarly, the experimental and calculated values of isentropic 
compressibility along with the percentage deviations at different mole 
fractions for toluene + o-xylene and toluene + aniline are reported in tables 
6.5 and 6.6 respectively. 

A close perusal of tables 6.1, 6.2, 6.3, 6.4, 6.5 and 6.6 reveals 
that the calculated values of u, p T and p s for both the mixtures show the 
same trend as experimental one. A close perusal of tables 6.1 and 6.2 
shows that ultrasonic velocity for both the mixtures increases as pressure 
increases at each mole fraction, and at a constant pressure, ultrasonic 
velocity increases as the concentration of toluene in the mixtures 
decreases. For both the mixtures, percentage deviation increases as 
pressure increases with maximum deviation of seventeen percent. 

From tables 6.3 and 6.4, it is evident that computed values of 
isothermal compressibility for both the mixtures decrease as the pressure 
increases, and at a constant pressure it decreases as the mole fraction of 
toluene in the mixtures decreases. The percentage deviation increases as 
the pressure increases, and maximum percentage deviation for both the 
mixtures is nearly fifteen percent. 
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A close observation of tables 6.5 and 6.6 reveals that the 
variation of isentropic compressibility with pressure and also with mole 
fraction shows the same trend as isothermal compressibility i.e. it 
decreases as pressure increases or as the concentration of toluene in the 
mixtures decreases. In this case also, the percentage deviation increases 
as pressure increases but it is greater than that as observed in the case of 
isothermal compressibility. Here, the maximum percentage deviation is 
nearly forty six percent for toluene + o-xylene and is nearly forty three 
percent for toluene + aniline. 



Table 6.2 



Table 6.3 


67 



Table 6.4 



Table 6.5 



Table 6.5 
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Introduction: 

The characteristic, both of medium and the ultrasonic wave 
propagating through it, is not affected by one another provided the 
propagating wave has infinitely small amplitude. But, the measurements of 
the ultrasonic propagation constants e. g. absorption and velocity become 
unreliable* 1 * when sound waves of finite amplitude propagate through 
fluids, because several interesting non-linear effects such as acoustic 
streaming, wave distortion, harmonic generation and sound saturation 
occur due to unequal attenuation suffered by its different components (the 
absorption coefficient is proportional to the square of frequency). 

These non-linear effects are, mathematically, represented by 
the non-linear terms in the equation of wave propagation. The non-linear 
departure of the propagation process is symbolized by the factor (B/A), 
called the non-linearity parameter. If the isentropic equation of state of a 
medium is represented as a Taylor series and only the quadratic and 
linear terms are considered, the ratio of the respective coefficients of these 
terms, B/A, is used as a parameter to describe the non-linearity of the 
medium. B/A of liquids can be obtained from the distortion of finite 
amplitude waves and the variation of sound velocity with temperature and 
pressure. 

The non-linearity parameter, B/A, plays a significant role in 
non-linear acoustics and its determination is of increasing interest in a 
number of areas ranging from under water acoustics to medical science. 
From the knowledge of the parameter of non-linearity, one can gain 
information about some physical properties (2,44,45) of the liquids such as 
internal pressure, intermolecular spacing, acoustic scattering and 
structural behaviour etc. 

A number of experimental and theoretical studies have been 
performed on the non-linearity parameter of liquids making use of 
phenomenological (2_4) or thermodynamic methods (5, 6) . Also, some studies 
on B/A have been performed from the view point of theory of liquids (7) . 
The values of B/A for a series of liquids have long been determined by 
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many authors using different methods (8 " 10) i but measurement of this 
parameter for biological media (11 ' 13) including living tissues is one of the 
more recent interest as it can provide important information for the 
ultrasonic application in biological research related to diagnosis and 
therapy. It was discovered that the non-linearity parameters of biological 
media differed strikingly from each other according to their molecular 
structure, so we can think that B/A will give some information of the 
molecular structure of a medium. But we know very little about the relation 
between them. 

Beyer et al (1,2,8,14) were the first who presented the theoretical 
formalism of non-linearity parameter and determined it for a number of 
organic liquids using sound velocity and other thermodynamic data. 
Hartmann et al (15) discussed the phenomenon (non-linear behaviour of 
liquids) in view of the potential energy between the molecules by 
assuming that sound velocity depended on volume only. Endo (5,7,16, -' ) 
obtained the expression for the non-linearity parameter in the form of a 
polynomial in terms of the ratio of specific heats and showed that it can be 
computed from thermodynamic coefficient of sound velocity. Hartmann 
and Balizer (15) have determined the values of non-linearity parameter for 
n-alkanes using a new equation of state. Theoretical evaluation of these 
non-linearity parameter for various class of liquids has been done by 
several workers (18 ' 20) Beyer (1 ' 2,8,14) , Nomoto (21) , Kor et al (22 ' 24) , Swamy 
et al (18) , Pandey et a | (25 ' 27,39) and Sharma (28 ' 32> succeeded in computing 
non-linearity parameter for a number of liquids using thermodynamic 
approach. Tong et al (33) proposed a simplified method for calculating the 
B/A values for pure organic liquids making use of Schaaffs equation (34) for 
sound velocity. Recently, using this method, Jugan et al (35) and Sanguri et 
al (36) calculated B/A values of binary liquid mixtures, and higher alkanes 
respectively. : • 

In the present work, computation of B/A for pure liquids has 
been carried out at elevated pressures by four different methods. 
Furthermore, a comparative study of B/A values obtained from these 
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methods has also been made in order to review above mentioned 
approaches, and the merits and demerits of these different methods are 
discussed in the light of molecular structure and intermolecular 
interactions. 

Thermal expansion coefficient (a) and isothermal 
compressibility (p T ) of pure liquids are evaluated from recently 143,46,1 
proposed relations, : 

75.6 *10 -3 

a ml 19 1/2 1/3 , 


17.1 * 10 -4 
Pt T A,9 p 4 % 2 


--(ii) 


Adiabatic compressibility (p s ) is calculated from the thermodynamic 
relation, 


&=-f- 

u P 


— (iii) 


The values of ultrasonic velocity and density are taken from the 
literature (37) . Reduced volume is obtained using eq. (2.1 1). 


Theory: 

Ballou (38) proposed an empirical rule for the non-linearity 
parameter using data from Beyer (8) on organic liquids, water, liquid metals 
etc. According to the empirical observation made by Ballou, there is a 
linear relation between the non-linearity parameter of liquids and 
reciprocal sound speed as, 


£- 0 . a + iLm 4 > 

A u 


--(7.1) 


where u is the sound velocity expressed in m/sec. 

Hartmann (38) has shown theoretically the physical basis for the 
above mentioned empirically observed relation between the non-linearity 
parameter of liquids and reciprocal sound speed assuming that the inter- 
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molecular potential energy is the dominant factor in determining sound 
speed, and its derivatives in liquids. This assumption implies that the 
thermal pressure contribution to the sound speed is negligible in 
comparison to the internal pressure contribution. The purpose of the 
theoretical calculation is to provide physical insight into the proportionality 
of B/A with 1/u. The expression for non-linearity parameter, B/A, due to 
Hartmann (38) is 

B=2 + (0 1 9SM0 4 ) _ (7 . 2) 

A u 

This theoretical result was obtained without using any adjustable 
constants and is in very good agreement with the empirical result. 

The non-linear behaviour due to the propagation of finite 
amplitude wave in fluids can be described by the equation of state P= 
P (p, S), where P is the pressure,# pis the density and S is the entropy of 
liquid. 

Non-linearity parameter, B/A, can be calculated from the 
modified thermodynamic equation of state (14) for liquids, namely, 


( 7 . 3 ) 

-(7.4) 

in which subscript ‘o’ indicates the state without perturbation. 

Since “» = fc)j p-p. . 

we get, 


P = P« + A 


P-P o 


, + A-Po 
Po J 2 l Po 


. (dP 

where 


and & ~ Po 


S p - po 
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and, A = p 0 uq 
therefore, 

£ _ Po [ Suq 


A u^} P = P0 


P o [ du Q ap 

ap 'dp 


- Po 


Po 


--(7.6) 


In the following deduction subscript ‘o’ will be omitted but all the 
variables will take the value at equilibrium. 

Thus, we can write from eq. (7.6), 


B 

7 =p 


—(7.7) 


here, it is clear that B/A is related to sound velocity. 

Tong et al (33) applied Schaaffs equation {34) for sound velocity 
in eq. (7) and obtained an equation for B/A as, 


-= J(o) + J(x') 
A 


where, ■/(„> = (i - l) i 

on simplification, 

.7-1 


J (°) = 
and J(x') = 


ccT 

2(3 - 2 x ') 2 


3 (jc' - 1)(6 - Sjc') 


—(7.8) 


—(7.9) 

—(7.10) 


here, y is the ratio of specific heats, p the density, a the thermal 
expansivity, (3 T the isothermal compressibility and x’ is the real volume of a 
molecule which is the ratio of molar volume, V, to b, where b is the van der 
Waals constant given by (34) 
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r i 

1 M 


\ 3/RT 

p) 


where M is the molecular weight, R the universal gas constant, u the 
ultrasonic velocity and P is the pressure. 

Because the liquid volume is always larger than the real 
volume of molecules, x’ is always larger than 1. The expression (33) for 
sound velocity contains a term (6-5x’) in the numerator. Since the 
ultrasonic velocity can never be negative, the value of x’ can never be 
more than 1 .2. Thus, we have the condition for x’ as 

1<x’ <1.2 -—(7.12) 

From eq. (7.10), it is evident that as the value of x’ tends to 1.2 the 
value of J(x’) tends to infinity which is an absurd result. Hence, value of x’ 
plays a very important role in determining B/A values using Tong-Dong 
method. 

General expression for the non-linearity parameter in terms of 
the acoustical parameters of liquids has been derived using the 
expression for the sound velocity (u), and introducing the contribution due 
to isobaric acoustic parameters (K) and the isothermal acoustic parameter 
(K’). The properties have been deduced with the help of the thermal 
expansion coefficient (a), whose dependence on temperature and 
pressure may conveniently be utilized to evaluate them. 

Moelwyn-Hughes parameter (40 ' 41) (C^is an important 
thermodynamic quantity in studying the thermodynamics of molecular 
surface and thermo-acoustical properties. It can be expressed as 



The Sharma parameter (S 0 ) may be expressed as (30) 
s ° ) 3+4ar) 
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= IS* (\ + 2aT)l K t 

= 0 + 2«rX3 + 4ar)/^ r --(7.14) 

where X is the isochoric temperature coefficient of internal pressure, 
r 2 (d In a) ] 


X = 


d In P, 


H 


(d In T) v 


d\nT )y 
-2(\ + 2aT / (V) c ' 



-(7.15) 


k T is the reduced isothermal bulk modulus, i. e. 

K r =(k t /k\ )=(V) C ' —(7.16) 


v is the reduced specific volume and expressed as, 
V(V/V*) = (S*/l+aT) 3 

{fyM-T 


where S =1 + 


m 


— (7.17) 
—(7.18) 


V* and k' t are respectively the hard core specific volume and bulk 
modulus at absolute zero temperature. 

Huggins parameter (F) may be expressed in terms of 
temperature and pressure derivatives of compressibility (26) , 


fdlnX T V 

CaT ' 

( d\n K t ^ 

[ d In T ] P 

+ UrJ 

[ d p X 


2aT 


f d\n K t \ 
\ d In T J v 


--(7.19) 
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using eqs. (7.14) and (7.19), the Huggins parameter may be 
expressed in terms of S 0 as, 


F + S* = 2^1 + (s*S 2 )/ir7’ j = 2^1 + ~T 


-(7.20) 


The isothermal acoustical parameter (K’) is related with the available 
molar volume (V a ) of substance by the expression, 


f = (K + lj =(K 


— (7.21) 


In this equation the isothermal (K), isochoric (K) and isobaric (K) 
acoustical parameters are given by, 


'HIT), 

_ ( d In u \ 
l d In V ) r 




(dXnKi 


\ d\nT J v 


2aT 


l 2 aT 


and k ■■ 


1 ¥ d In i< ^ 
a ){ dT ) P 


= i[i+^(i + ar)/ar] 


-(7.22) 


--(7.23) 


—(7.24) 


Beyer’s parameter of non-linearity, (B/A), which is a particular 
combination of the temperature and pressure derivatives of the sound 
velocity may be expressed as (42) , 

f 2Mu 2 aT V d In n \ 


B _ 

'2 Mu 2 ) 

( din u V 

A ~ 

< v J 

l d P )> 


■T 


/ Vlna T 

t dT % 


--(7.25) 



= 2 yK' - 2K(y - l) 

= 2K + 2yK" — (7.26) 

where M is the molecular weight, y = C P /C V the heat capacities ratio, 
and C P and C v are respectively the isobaric and isochoric heat capacities 
of the substance. 

Results and Discussion: 

The values of non-linearity parameter, B/A, of pure liquids 
have been evaluated at elevated pressures and at temperature of 303.15 
K with the help of four different methods: Ballou’s empirical relation 
(eq.7.1), Hartmann’s theoretically derived relation (eq. 7.2), Tong - 
Dong method using Schaaffs equation for sound velocity (eq. 7.8) and 
Sharma's method vide eq. (7.26). The three pure liquids under the present 
consideration are: toluene, o-xylene and aniline. The different parameters 
of these liquids, needed for the computation, have been taken from table 
3.1. Experimental values of ultrasonic velocity and density are taken from 
the literature (37) . 

The calculated values of non-linearity parameter, B/A, for 
toluene, o-xylene and aniline, using the above mentioned methods, are 
reported in tables 7.1, 7.2, and 7.3 respectively. The parameters x', J(o), 
J(x’) for toluene, o-xylene and aniline used in eq. (7.8) are also reported in 
the corresponding tables. A close perusal of tables 7.1, 7.2 and 7.3 
reveals that the values of B/A calculated using eqs. (7.1) and (7,2) show a 
decreasing trend with increase of pressure, whereas the B/A values 
obtained from eq. (7.26) increase with increasing pressure. However, the 
variation in the B/A values due to this method is very small. Table 7.1 
shows that B/A values for toluene, using eq. (7.8), generally, decrease: 
with increase of pressure. In this case, some values of B/A are found to be 
negative. It is because of the value of x’ which is greater than 1.2, the 
maximum limit of x’ in the Tong-Dong method. In fact, in the Tong-Dong 
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mehtod, the B/A values are highly sensitive to the values of x’. It has been 
noticed that even small changes in values of x’ give very absurd values of 
B/A. This is mainly due to the factors (3-2x’) in the numerator and (6-5x') 
in the denominator in the mathematical expression for J(x’). in fact, good 
values of J(x’) or B/A result when x’ lies between 1.10 and 1.14. It is clear 
from table 7.2 that B/A values for o-xylene decrease with increase of 
pressure. Table 7.3 reveals that, in the case of aniline, there is a decrease 
in B/A values with increase of pressure upto 70 MPa, but after this value of 
pressure, B/A values show an increasing trend with pressure. 
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Introduction: 

Study of non-linearity parameter, B/A, of liquid mixtures {1 ' 9) is 
important in view of the information it may yield on the interaction in liquid 
mixtures. B/A values of the liquids have been interpreted as the quantity 
representing the magnitude of the hardness of the liquid (10) which may be 
considered to be true for the liquid mixtures as well. 

Sehgal and Porter (11) determined the S/A values for alcohol- 
water mixtures and Endo (12) calculated B/A for benzyl alcohol-ethanol 
system. Jugan (13) used the Tong -Dong method to calculate B/A for a 
number of binary liquid mixtures. Recently, Kuchhal et al (14) and Pandey et 
al (15,16) determined B/A values for pure liquids and binary liquid mixtures. 

Although detailed studies of the non-linearity parameter of 
binary liquid mixtures have been reported, to the best of our knowledge, 
computation of B/A for liquid mixtures at elevated pressures has not been 
done so far. 

In the present chapter, values of non-linearity parameter for 
binary liquid mixtures have been computed at elevated pressures applying 
four different methods (9,17 ' 19) , and a comparative study is made. The 
thermodynamic parameters required to calculate B/A values of the liquid 
mixtures, are obtained from the Flory’s statistical theory as discussed in 
chapter-2. Experimental values of sound velocity, density and other 
related properties of two binary mixtures (toluene +i o-xylene & toluene + 
aniline) under taken for the present investigation, have been taken from 
the paper of T akagi et al {20) . 

Theory: 

Equations (7.1), (7.2), (7.8) and (7.26) have been used to 
calculate the values of non-linearity parameter, B/A, of binary liquid 
mixtures under the present consideration at elevated pressures. The 
theoretical method employed has already been detailed out in chapter 


seven. 


Results and Discussion: 


The two binary liquid mixtures under the present investigation 
are: toluene + o-xylene and toluene + aniline. The non-linearity parameter, 
B/A, of these liquid mixtures are calculated at elevated pressures ranging 
from 0.1 MPa to 160 MPa and at temperature of 303.15 K for four different 
mole fractions. Experimental values of ultrasonic velocity and density of 
the liquid mixtures are taken from the literature (20) . 

For the parameters of liquid mixtures, needed for present calculation, we 
have used the Flory’s theory in a manner as described in chapter-2. 

The calculated values of B/A for toluene + o-xylene using 
different methods are reported in table 8.1. The values of x’, J(o) and J(x') 
used in Tong-Dong method at different mole fractions for this system are 
reported in table 8.1(a). Similarly, table 8.2 enlists the values of B/A 
calculated through different methods for toluene + aniline system. Table 
8.2(a) contains the values of x’, J(o), J(x’) for toluene + aniline. 

A close perusal of tables-8.1 and 8.2 reveals that the B/A 
values calculated by Ballou’s empirical relation and Hartmann’s 
expression decrease as the pressure increases but the values of B/A 
computed by Sharma method show an opposite trend i. e. it increases as 
the pressure is increased. B/A values calculated by Tong-Dong method 
initially decrease with pressure upto a certain pressure, and after that it 
increase with the increase of pressure for both the mixtures. The value of 
pressure at which the trend changes decreases as the composition of 
toluene decreases in both the mixtures. Also, the non-linearity parameter 
decreases, in general, with the decrease in concentration of toluene in the 
mixtures upto a certain pressure. After this pressure, as the mole fraction 
of toluene in the mixtures decreases, the B/A values tend to increase. 
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